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Abstract. The lithium abundances in metal-poor halo stars are of importance for cosmology, galaxy evolution 
and stellar structure. In an attempt to study possible systematic errors in the derived Li abundances, the line 
formation of Li I lines has been investigated by means of realistic 3D hydrodynamical model atmospheres of halo 
stars and 3D non-LTE radiative transfer calculations. These are the first detailed 3D non-LTE computations 
reported employing a multi-level atomic model showing that such problems are now computationally tractable. 
The detailed computations reveal that the Li I population has a strong influence from the radiation field rather 
than the local gas temperature, indicating that the low derived Li abundances found by Asplund et al. (1999) 
are an artifact of their assumption of LTE. Relative to 3D LTE, the detailed calculations show pronounced over- 
ionization. In terms of abundances the 3D non-LTE values are within 0.05 dex of the ID non-LTE results for 
the particular cases of HD 140283 and HD 84937, which is a consequence of the dominance of the radiation in 
determining the population density of Li I. Although 3D non-LTE can be expected to give results rather close 
(~ ±0.1 dex) to ID non-LTE for this reason, there may be systematic trends with metallicity and effective 
temperature. 
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1. Introduction 

Lithium abundances of mctal-poor halo stars have had a 
prominent place in astrophysics and cosmology ever since 
the first discovery of Li in such stars (Spite & Spite 1982) 
due to their wide-ranging implications. Firstly, the uni- 
formity of the Li abundances in halo stars (Ryan et al. 
1999) presumably reflects the primordial Li abundance 
stemming from Big Bang nucleosynthesis, thus allowing an 
estimate of the baryon density of the Universe. Secondly, 
any observed slope in the Spite-plateau with metallicity, 
when combined with the corresponding behaviour of ^Li, 
Be, B and O, can be interpreted in terms of Galactic chem- 
ical evolution, in particular the history of cosmic ray spal- 
lation in the early Galaxy. Thirdly, since it is a fragile 
element, Li can function as a tracer of stellar mixing; the 
thinness of the Spite-plateau certainly severely limits the 
allowed amount of Li-depletion in halo stars (Ryan et al. 
1999). 

In order to extract Li abundances from an observed 
spectrum it is necessary to have suitable models of the stel- 
lar photosphere and the line formation process. Since tra- 
ditional abundance analyses utilize ID hydrostatic model 
atmospheres based on local thermodynamic equilibrium 
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(LTE) and a rudimentary description of convection in the 
form of the mixing length theory (MLT), there may be sig- 
nificant systematic errors, which could distort the derived 
conclusions. The detailed line formation, including depar- 
tures from LTE but in ID hydrostatic model atmospheres, 
were investigated by Carlsson et al. (1994) who found rel- 
atively small non-LTE abundance corrections (^ 0.05 dex) 
for low metallicity stars. To investigate the effects of con- 
vection on the line formation, Asplund et al. (1999) ap- 
plied the new generation of 3D, hydrodynamical model 
atmospheres to abundance analyses of metal-poor stars 
assuming LTE. Due to the much lower temperatures en- 
countered in the line-forming layers than in classical ID 
model atmospheres, the derived abundances differ drasti- 
cally. Most notably, it was found that the primordial Li 
abundance may previously have been over-estimated by 
~ 0.3 dex. Although such 3D model atmospheres should 
be a more realistic description of the stellar photospheres 
(Asplund et al. 2000), the simplification of LTE in the line 
formation is a crucial assumption. Asplund et al. (1999) 
indeed warned that the steep temperature gradients may 
be prone to significant over-ionization of species like Li i. 
The aim of the present paper is to investigate possible de- 
partures from LTE by performing realistic 3D non-LTE 
calculations for Lii lines, while we defer a detailed in- 
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vestigation of the line formation to a subsequent article. 
Preliminary calculations have been presented in Asplund 
(2000) using a smaller Li atom and without consideration 
of line-blanketing. 

2. 3D non-LTE radiative transfer of Li I lines 

From the ab-initio 3D radiative-hydrodynamical convec- 
tion simulations of the metal-poor halo stars HD 140283 
and HD 84937 presented in Asplund et al. (1999), two 
representative and independent snapshots have been se- 
lected from each simulation for the non-LTE calculations. 
Additionally, two snapshots were taken from a similar so- 
lar simulation (Asplund et al. 2000a) for comparison pur- 
poses. The original simulation data cubes with dimensions 
100^x82, were interpolated prior to the non-LTE calcu- 
lations to a 25^ X 100 grid with an improved vertical reso- 
lution. Test calculations with 50^ x 100 snapshots verified 
that the procedure had insignificant effect on the 3D non- 
LTE abundance corrections. Typically the atmospheric 
structures cover the region —6 < logrsoo < 2. Further 
details of the numerical simulations are available in Stein 
& Nordlund (1998) and Asplund et al. (1999, 2000). 

The Li model atom employed in the present study 
is identical to the one compiled by Carlsson et al. 
(1994). The adopted 21-level atom consists of in total 
70 bound-bound and 20 bound-free radiative transitions. 
We feel confident that the employed Li atom is suffi- 
ciently extended as tests with a more restricted 12-level 
atom revealed non-LTE abundance corrections only about 
0.03 dex larger than with the 21-level atom for HD 140283 
and HD 84937. The input data for level energies, radiative 
and electron coUisional transitions are in all cases the same 
as in Carlsson et al. (1994); accordingly, collisions with 
neutral hydrogen have not been included. We have verified 
that even inclusion of H collisions according to the classi- 
cal Drawin (1968) recipe has a very minor impact on our 
results; for HD 140283 the abundance corrections are de- 
creased by 0.05 dex and 0.01 dex for multiplication factors 
of 1.0 and 0.01, respectively, of the standard Drawin esti- 
mates. Additional atomic data such as background opaci- 
ties (including line-blanketing for photo-ionization transi- 
tions) are taken from the MARCS package (Gustafsson et 
al. 1975 and subsequent updates). 

The non-LTE calculations for Li have been performed 
with multi3d (Botnen 1997; Botnen & Carlsson 1999), 
which is essentially a 3D- version of the widely used multi- 
code for ID statistical equilibrium problems (Carlsson 
1986). multi3d iteratively solves the rate equations with 
a consistent radiation field obtained from a simultane- 
ous solution of the radiative transfer equation along 24- 
48 inclined rays. The formal solution of the radiative 
transfer equation is computed using a short character- 
istic technique, making use of the horizontal periodic 
boundary conditions of the hydrodynamical model atmo- 
spheres. The Doppler shifts introduced by the convective 
motions are taken into account and therefore no micro- 
and macroturbulence enter the analysis (Asplund et al. 
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Fig. 1. Grotrian termdiagram for the adopted 21-level Li 
model atom. All levels are connected with the Li ii ground 
state by photo-ionization transitions. 

Table 1. Predicted LTE and non-LTE flux equivalent 
widths (in pm= 10 mA) of the Lii 670.8 nm line for the 
Sun, HD 140283 and HD 84937 in ID and 3D. The 3D 
equivalent widths are the average of two independent 
snapshots. 
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^NLTE 


Sun 


1.1 


0.40 


0.34 


0.55 


0.37 


HD 140283 


1.8 


2.38 


2.18 


3.84 


1.96 




2.2 


4.90 


4.76 


7.67 


4.21 


HD 84937 


2.0 


1.31 


1.44 


1.79 


1.11 




2.4 


2.91 


3.29 


4.00 


2.55 



2000). Typically ^ 6 accelerated lambda-iterations with 
linearization and preconditioning of the rate equations 
(Scharmer & Carlsson 1985) were necessary to achieve the 
stipulated convergence criteria max(i5ni/ni) < 10"''. The 
resulting angle-dependent radiation field was subsequently 
used to compute flux profiles. Various test calculations 
with 3D-extended plane-parallel homogeneous model at- 
mospheres ensured that the same results were obtained 
with multi3d and multi. 



3. Departures from LTE in Li I line formation 

Table ^ summarizes the resulting 3D LTE and non-LTE 
line strengths for the different snapshots and input Li 
abundances together with the corresponding ID results. 
Clearly, there are very significant departures from LTE, 
in particular for the two metal-poor stars. As predicted 
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Fig. 2. Granulation pattern in snapshot 1 of HD 140283 
seen in disk-center continuum intensity (upper left panel) 
and equivalent width of the Lii 670.8 nm line in LTE [up- 
per right panel) and non-LTE [lower left panel)] the equiv- 
alent width images have the same relative intensity scale 
to emphasize the overall difference in line strengths. Also 
shown is the ratio of the non-LTE and LTE equivalent 
widths [lower right panel). 

by Asplund et al. (1999), over-ionization plays the domi- 
nant role. 

Fig. 12 shows images of the continumn intensity and the 
equivalent width of the Li i 670.8 nm line in LTE and non- 
LTE at disk-center (/^ = 1.0). The characteristic granula- 
tion pattern with warm (bright) upflows and cool (dark) 
downflows is clearly visible in the continuum intensities. 
In LTE the Li line typically is strong above the granules, 
a consequence of the low temperatures in the optically 
thin layers due to the temperature contrast reversal in 
the convectively stable layers, which translates to a high 
population density of Li i. Over the downflows the temper- 
atures tend to be higher than average, thus making the Li 
line weaker. In non-LTE quite the opposite happens: in 
general the line is weaker above the upflows, which im- 
plies that the population density is largely controlled by 
the radiation field rather than the local gas temperature. 

The same conclusion is obvious from Fig. 13 which 
quantifies the behaviour in Fig. |3 The Li i lines tend to 
increase in strength with continuum intensity in LTE but 
with a very large scatter, since the temperatures in the 
higher atmospheric layers are not one-to-one correlated 
with the temperatures in the continuum-forming layers 
below. In non-LTE the scatter is much smaller for regions 
with high continuum intensities, which again emphasizes 
that the line formation process is largely governed by the 
non-local properties of the radiation field. The 670.8 nm 
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Fig. 3. Predicted intensity (^ = 1,0) equivalent widths of 
the Li I 670,8 nm (left panel) and 610,4 nm (right panel) 
lines across the granulation pattern in LTE (upper panel) 
and non-LTE (middle panel) for the HD 140283 simulation 
with logeLi = 2,20. The mean intensity equivalent widths 
are denoted by horizontal lines. Also shown are the ratios 
of non-LTE and LTE equivalent widths (lower panel). 



line is, however, weakest in relatively low continuum inten- 
sity regions (/™nt/^jcont^ ^ g.g in Fig, EJ) which are im- 
mediately adjacent to very bright regions: the non- vertical 
hot radiation field from neighboring granules cause signifi- 
cant over-ionization (low line opacity) which when coupled 
to the shallow temperature gradient of downflows makes 
the line very weak (e,g, at a; = 21 Mm and y = 18 Mm in 
Fig, EJ, The Lii line typically is strongest in the mid- 
dle of larger areas of downflowing material, which has 
weaker photo-ionizing radiation field (e,g, at x = 9 Mm 
and y = 13 Mm), The non-LTE behaviour of the Lii 
670.8 nm line shown in Figs. |21 and |31 is similar to the 
findings of Kiselman (1997, 1998) and Uitenbroek (1998) 
in the case of the Sun, which are confirmed by spatially 
resolved solar observations (Kiselman & Asplund 2001). 
The Lii 610.4 nm line show a similar behaviour as the res- 
onance line although less dramatically due to the weakness 
of the line and its higher excitation potential (Fig.lSJ. 
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Fig. 4. Variation of the departure coefficient of the Li i 
ground level Ni/N^^^ (upper panel) and S^^-^/B^T) 
(lower panel) as a function of atmospheric height for the 
HD 140283 simulation. Selected vertical columns arc con- 
nected with solid lines. 



Further support for the departures from LTE being es- 
sentially an over-ionization effect comes from the variation 
with height of the departure coefficient of the Li i ground 
level (Fig. 0J). From the LTE-expectation in the deeper 
layers, the population density falls rapidly with height as 
Li I is photo-ionized by radiation from below. The non- 
LTE line formation of the Lii 670.8 nm line is mainly an 
optical depth effect since the line source function S"^™-^ 
in general remains reasonably close to the Planck func- 
tion By{T) (Fig.^. As in the ID case, the over-ionization 
is predominantly driven by photo-ionization from the 2p 
state for which Jy/Bi, is near its maximum (Carlsson 
et al. 1994). The inclusion of line opacities to the back- 
ground opacities for bound-free transitions is thus of im- 
portance: without line-blanketing the non-LTE abundance 
corrections for HD 140283 would be about 0.07 dex higher. 
As expected the radiative rates dominate over the corre- 
sponding coUisional rates for the 2s and 2p levels. 

The speculation by Kurucz (1995) that the strength of 
the Li I 670.8 nm line could be seriously overestimated with 
ID models due to substantial departures from LTE in the 
presence of temperature inhomogeneities, is thus partly 
validated by our results, although for the wrong reasons 
and to a much lesser extent. Due to his simple two-stream 
MLT-convection approach, Kurucz also failed to recognize 
the compensating effect due to convection (Asplund et 
al. 1999; Asplund 2000). Cayrel & Steffen (2000) reached 
similar conclusions to ours based on calculations using a 
5-level Li atom and a 2D simulation of a metal-poor Sun. 



4. Discussion 

With the aid of Table it is possible to interpolate the 
computed ID and 3D line strengths to the observed equiv- 
alent widths of Li I 670.8 nm for HD 140283 (4.65 pm) 
and HD 84937 (2.44pm). It is noteworthy that the 3D 
non-LTE abundances only differ by ^ 0.05 dex from the 
corresponding ID non-LTE results for HD 140283 and 
HD 84937. The 3D LTE abundances on the other hand 
are distinctly different by 0.2 — 0.3 dex. As cautioned by 
Asplund et al. (1999), LTE is clearly a very poor assump- 
tion for Li I in 3D models of halo stars with their low 
surface temperatures and steep temperature gradients ^ . 

The 3D non-LTE abundances are very similar to the 
ID non-LTE results. The main reason for this is that the 
line strength is determined by the opacity which is set 
by the photoionizing radiation field while collisions are 
found to be relatively unimportant according to our cal- 
culations^. This radiation is optically thin at the height 
of formation of the Lii resonance line and is thus equal 
to the observable emergent radiation field. Since ID and 
3D models produce similar (±10%) emergent intensities at 
the relevant wavelengths (Asplund & Garcia Perez 2001), 
the photo-ionization radiation fields will not differ signif- 
icantly. The final line strengths will however also depend 
on the details of the temperature structures in the deep 
atmospheric layers which differs between ID and 3D mod- 
els. 

From the results obtained here for two typical halo 
stars we expect that the 3D non-LTE abundances will be 
similar (« ±0.1 dex) to the ID non-LTE results (and to ID 
LTE given the small ID non-LTE abundance corrections 
according to Carlsson et al. 1994) in general for halo stars. 
However, given the great interest in accurately determin- 
ing the primordial Li abundance, the intrinsic scatter in 
the Spite-plateau and the possible existence of a trend 
with metallicity for the Li abundances (e.g. Ryan et al. 
1999), it is clearly important to extend the 3D non-LTE 
calculations to additional stars with different parameters 
to determine exactly the net 3D effects. We are currently 
working towards this goal (Asplund et al., in preparation). 
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